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The ground state of the multi-flavor 2-dimensional quantum electrodynamics (QED2) is deter- 
mined in the presence of finite baryon density, and it is shown that the model possesses two phases, 
the low density phase where the external baryon density is totally screened, and the high density 
phase, where the screening is only partial. The renormalization of the bosonized version of the 
model is also performed for both the zero and the finite density model giving massless multi-flavor 
QED2 in both cases. 
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O ' I- INTRODUCTION 

o : 

The phase structure and the confinement mechanism of non-abelian models are usually investigated in much simpler 
,-0 ■ usually two dimensional 'toy' models [l|,|2|. These models are sometimes analytically solvable, e.g. the massless QED2 
rS^ I (usually referred to as Schwinger model) shows the chiral condensate and a mass gap which are supposed to be 
essential elements in modern physics. The bosonized version of the Schwinger model is a free theory, which turns into 
' interacting when the fermionic mass is nonzero QED2 possesses a phase transition at rn/g c ~ 0.31 as was shown 
by density matrix renormalization group (RG) technique [5| or by continous RG method [6(. The critical value of 
m/g c separates the large coupling (g 3> m) phase with a unique vacuum characterized by the field variable = 0, and 
the weak coupling (g -C m) phase where due to spontaneously broken reflection symmetry the model has non-trivial 
O .■ vacua, at around <f> = ±y/ir/2. 

The ground state of quantum chromodynamics (QCD) is a color superconductor at high densities Q and a periodic 
chiral condensate appears in the coordinate space at large number of colors N c Q which was recapitulated in the 
framework of renormalizaton group (RG) equations by mapping QCD into a Thirring type model Q. The high 
, density behavior of non-abelian models are investigated in the framework of toy models too [icil ]. The Schwinger 
' model remains exactly solvable in the presence of an external finite charge density and also shows periodic chiral 
&\ ', condensate [HI, EH • -"- n our previous work jj| we showed that in the finite density QED2 the system exhibits a single 
periodic phase in the thermodynamic limit for arbitrary charge densities, furthermore the periodic structure built up 
in the ground state has decreasing amplitude and wavelength with increasing charge density. It is assumed 2] that 
a phase transition appears in the 2-flavor QED2 as the density is increased. By the bosonization technique [14] both 
the massless and the multi-flavor QED2 can be converted to such local scalar field theoretical models which contain 
periodic self-interaction potentials. The bose form of the multi-flavor QED2 is a sine-Gordon-type scalar field theory 
[Tol liH [ijj , where periodic self-interactions are described by 2-dimensional sine-Gordon (SG) fields which coupled 
by a mass matrix. This coupled SG (or layered sine-Gordon (LSG)) model can also be used to describe the vortex 
dynamics of magnetically coupled layered superconductors [l8j , where the number of flavors in the high energy model 
is identical to the number of layers of the condensed matter system [l7j . Moreover coupled SG type models have been 
■ used to investigate the vortex dynamics of Josephson coupled superconductors [l9| . 

The phase transition of these models was obtained from the microscopic theory which is formulated in the high 
energy or ultraviolet (UV) region. In order to obtain the low energy/infrared (IR) physics, where the measurements 
are performed and the quantum fluctuations are taken into account we need renormalization. The original, fermionic 
models so as the toy models contain strong couplings and it disables performing a perturbative renormalization, and 
it makes rather difficult to develop a functional renormalization group (RG) method in the fermionic models since 
the evolution should be started from a perturbative re gion where the theory is almost interaction free. Furthermore 
the RG equations has to preserve the gauge symmetry [20j • However the bosonized version of the toy models can be 
easily treated by the functional RG method [g, Ht], HJ • 

Our goal in this article to clarify the phase transition of the finite density 2-flavor QED2 by the methods used 
in We calculate numerically (and also analytically in the low and high density limits) the ground state field 

configurations for finite external baryon density in the tree-level. The results show that an induced baryon density 
appears which screens the external one totally (partially) for low (high) external baryon densities respectively, showing 
the existence of two phases as is conjectured in Q. We also perform the RG procedure to determine the IR physics of 
the zero and finite density 2-flavor QED2 and then generalize the results for arbitrary number of flavors. We choose 
the so-called Wegner-Houghton (WH) RG method [22| in order to obtain the blocked potential for the model, which 
uses a gliding sharp cut-off k. The external baryon density produces coordinate dependent non-trivial saddle points, 
therefore we use the tree- level blocking relation @, [2l|, [23[ to get the blocked interaction potential. The RG evolution 
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gives the flow of the couplings for the LSG model from which the evolution of the couplings of the fermionic model 
can be easily determined. 

The paper is organized as follows. In Sect.[n] we show the connection between the fermionic and the bosonized 
models and then in Sect. lIIII we determine the ground state field configuration of the finite density 2-flavor QED2 and 
map its phase structure. We derive the evolution of the coupling in the framework of the WH-RG method in Sect. lIVl 
and determine the flow of the couplings in the case zero, low and high densities. Finally, in Sect. |V] the conclusion is 
drawn up. 



II. THE MODEL 



The multi-flavor QED2 containing TV Dirac fields with identical fermionic charge e and mass m has the Lagrangian 
density 



^ JV N 



(1) 



i=l 



where Fm = 8qAi — d\Ao. One can turn from fermionic field variables tpi, ipi into bosonic ones <pj by the bosonization 
rules dli| 



-cmM cos(2v / 7T</)i) 
- cmM sin ( 2 \pxfa ) 

^=£aud"fa 
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where N m means normal ordering with respect to the fermion mass m, M = e/y/w, and c = exp(7)/27r, with the 
Euler constant 7 = 0.5774. The presence of a non-vanishing external or background densities does not affect these 
transformation rules fill ] . The Hamiltonian of the system in Coulomb gauge is given by 



A f - e 2 f 

Ti = } y / ipi (x) (vyidi + m)^i(x) - — / jo, x \x 



y\jo, y , 



(3) 



with J x = dx° J^ L dx 1 and 



JV JV 

jo,x =■ ^2i>i(x)ipi(x) := —^9x^2 



fa(x). 



(4) 



The resulting bosonized form of the Hamiltonian is 

JV , JV 



n = N n 



JV 1 JV 2 / N \ N 

2 E n * (*) + 2 E^^w) 2 + h ( ? Mx) ) ~ cm2 T, cos (v^(*)) 



i=l 
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(5) 



where H{(x) denotes the momentum variable canonically conjugated to (f>i(x). Let us assume that the two distinct 
flavored fermion has the same mass m but opposite charge e. The resulting bosonic Hamiltonian corresponds to the 
LSG model 



H = N M 
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Mfa 



fa) 2 - u(cos(2\/nfa) + C0s(2v / 7T'J>2)) 



(6) 



with two layers (N = 2), and u = cme/y/ir. The assumption for the charges corresponds to the situation where the 
net electric charge is zero which means that the matter is colorless. It is assumed that this kind of colorless matter 
might exist in nuclear stars. Let us denote the two types of bose field as fa with i = 1,2. After introducing the new 
fields as 
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the charge density jo = p and the baryon density B can be written as 



P : 



-.On 



B = 



-du 



(8) 



An external baryon charge density can be introduced by replacing p by p + po. If we define the corresponding classical 
potential as 



Po 



=dic 



(9) 



then the uniform, constant external baryon density is simply <f> c — bx. We separate this linear term in the space 
direction cj> + = <p + + bx, introduce the chemical potential p and assume that it is non-vanishing on the interval 
[— L;L], Taking b = — 2/i2/y / 7r then Hamiltonian reduces to 



H = N 



M 



+ \(di<j>-) 2 + M 2 (j) 2 _ + ^n 2 + + i(<9i0+) 2 - u' cos(V2^<j)-) cos(V2^(4>+ + bx)) 



due to the normal ordering with respect to the bosonic mass M where the dimensionful coupling v! 



(10) 



(11) 



is introduced, and d = c(2/n) 1 ^ 4 . However, the periodic part of the Hamiltonian also gives a contribution to the 
mass spectrum, i.e. the normal ordering can also be defined with respect to the 'total' mass. In this case the relation 
between the mass gap of the multi-flavor bosonic model and the parameters of the original fermionic theory is 



M gap = 2.008 • m 7r + T e^+ T 



(12) 



where N is the number of flavors [Ta|. The Hamiltonian (fT0|) is a functional of the field configurations 0_ and <f>+, 
and it gives the energy of the system on the tree level. We look for the ground state field configuration of the model 
as the function of the density b. 



III. THE TREE LEVEL PHASE STRUCTURE 



The minimum of the energy is searched numerically among the static field configurations, ir±(x) = by means of 
conjugate gradient method as a function of the finite baryon density b. We made numerical calculations in order to 
minimize the tree level energy as the functional of the field variables tf>- and </> + . The results show that at the energy 
minimum cj>_ (x) = for all value of b. This result is not surprising since in the Hamiltonian in Eq. (|10p one has a 
massive sine- Gordon (MSG) model [J BUI, El 

with zero density for the field variable <p^(x), and this model exhibits 
trivial field configurations on the tree level approximation 13]. Therefore we can treat a simpler Hamiltonian 

n = l \*\ + \( d ^+? - u ' cos(\/27(0+ + bx)) , (13) 
which corresponds to a SG model [2|| with finite density. Eq. (fl"3"|) leads to the equation of motion 

(Qo - df) 4>+ + VW sin(\/2^> + +bx))= 0. 



(14) 

Considering the static equation of motion and using a simple redefinition of the field variable v / 27r(</>+ + bx) — > <j>+ 
gives the static SG equation (2?| which is identical to the evolution equation of a pendulum. Depending on the initial 
energy C of the pendulum, the model has two phases. When C is large then it makes periodic rotation corresponding 
to the kink (or antikink) crystal in the original model [2^ |. In the low energy phase the pendulum swings, and we 
have a kink-antikink (KA) crystal. The external baryon charge contribute to the energy by the term irb 2 . Then the 
large density limit corresponds the KA crystal, while at low densities the kink crystal solution appears. The general 
analytic solution is 



-bx 



V2? 



2ir 



(15) 
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with am(i,r 2 ) the Jacobian elliptic function and r — \j2j (2 + C/2iru'). The coupling u' plays key role in the RG 
procedure, so we keep it explicitly. We consider the limiting low density (r — > oo or small values of b) and high density 
(r — > or high values of b) cases. The series representation of the function am(i, r 2 ) is 



am(x, r 2 ) 



TITTX 



2K{r) ^ n(q(r) 2n + 1) if(r) ' 



(16) 



where K(r) is the complete elliptic integral of the first kind and q(r) = exp(— irK(l — r)/K(r)). In the low density 
limit q{r) —> — 1 and the second term is just the Fourier expansion of the first term with opposite sign in the r.h.s. of 
Eq. ([Ill]) , giving am(x, r 2 ) — > 0, so 



-bx 



(17) 



as is conjectured in [2j. It gives a constant induced density equal to the external one with opposite sign, resulting 
total screening. At large densities the field configuration changes. When b is large C ~ -Kb 2 and the amplitude A n 
of the nth Fourier expansion in Eq. (fT6|) scales as A n ~ 1 jb 2n therefore it is a good approximation to keep only the 
fundamental mode. Then the field configuration becomes 



/2nb 2 



sin(v27r6x), 



(18) 



a sinusoidal type one, giving sinusoidal induced density which only partially screens the external density, and the net 
baryon density turns to finite values signalling a phase transition going from low densities to high ones. In Fig. [T] two 
typical numerically obtained field configurations are plotted. The wavelength £ is insensitive for the fermion mass 
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FIG. 1: The form of the field configuration for b — 0.1 (upper figure) and 6 = 5 (lower figure). 



m, furthermore it decreases according to a power law in b for several orders of the density. The numerical fit give 
that I = \[2~k jb equal to the 'specific volume' of the external density similarly to the case of MSG model [HI, as 
can be seen in Fig.[5J The amplitude A is a constant and independent of fermion mass for low densities since it is 
independent of the coupling u' . Due to the linear form of its value is A = M/2 = yjn /2. Going to high density 
regime the amplitude decreases as A ~ 1/b 2 and an m dependence appears. We plotted the amplitudes A n in Fig. [3] 
The limiting low and high density cases, is in a good agreement with its analytically calculated value in Eqs. (fTT)) 
and (fTB")) . Fig. [3] shows that the phase transition appears around the critical value b c . Using the field configurations 
obtained in Eqs. (| 1 T[) and (fTB")) one can calculate the total energy density £ by inserting them back into Eq. (|13p . For 
low densities it is (I refers to iinear field configurations) 

b 2 
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FIG. 2: The baryon density dependence of the wavelength I of the field configuration 




FIG. 3: The baryon density dependence of the Fourier modes A„, n = 1,2, 3. An follow a power law behavior at high densities 
as it was obtained analytically. 



while in high densities one obtains (p refers to periodic field configurations) 



u' 2 



which shows that the energy density is always negative and the field configuration in Eq. (jT5J) is energetically favourable 
in comparison with a trivial one which would give £ = 0. In Fig.|¥]we plotted £ p and Ei and the numerically determined 
total energy density. At low densities the linear field configuration is preferable and for high densities the sinusoidal 
one. The high density phase then contains a coordinate dependent ground state which is supposed to flatten out 
when one takes into account the quantum fluctuations beyond the tree level approximation ||. The critical density 
b c which separates the low and the high density phases can be read off correctly from Fig. [3] as the intersection of 
the Fourier modes, and it is b c « 0.3. The baryon density at the origin B(x=0) = Bo can play the role of the order 
parameter of the phase transition. It is zero in the low density phase and is finite in the high one, as can bee seen in 
the inset of Fig.UJ According to the inset of Fig. 2] the phase transition seems a continous one. 
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FIG. 4: The baryon density dependence of the energy density in case of linear Si (dotted line) and periodic £ p (dashed line) 
configurations. We also plotted the numerically determined total energy density by a solid line. At low densities the linear, at 
high densities the periodic field configuration coincides with the numerical results. The inset shows how the baryon density Bo 
increases as the background baryon density b grows. 



IV. RENORMALIZATION 



The phase structure of the finite density 2-flavor QED2 is determined at the tree level where the observational 
scale k is in the UV regime. We use the RG method to lower the scale k into the IR limit in order to consider 
the effect of the quantum fluctuations systematically. Here we perform the RG treatment of the LSG model for 
zero, and then for finite density. We treat the model in a more general form containing the arbitrary frequency 
parameter (3. The bosonization works for the specific choice 1 = 4tt and gives the Hamiltonian in Eq. ([6]), and in this 
article we map the phase structure as the function of the baryon density b. However the phase structure can be also 
considered in terms of the frequency (3 with the critical value /3 2 = 167r in the LSG model [2a |. At the scale where the 
microscopic theory is formulated the local potential approximation (LPA) is usually valid, implying that during the 
blocking transformation a constant field configuration gives the minimum of the potential [22L l29j , furthermore the 
RG methods based on the evolution of the effective action [30] are also formulated with constant field configurations. 
We showed in the preceding section, that for the finite-density 2-flavor QED2 it is not the case. The non-trivial 
saddle point field configuration makes unstable modes, and one should take into account coordinate dependent field 
configuration to find the real extremum of the Lagrangian, or at least better than the constant field configuration 
could give [23|, [3l[ . Therefore we choose the WH-RG method where one can change to tree level blocking relation 
when the non- trivial saddle point appears (23|. In this RG method the blocked potential can be derived by means 
of the differential RG in momentum space at the gliding momentum scale k £ [0, A]. In the WH-RG approach the 
field fluctuations are ordered by their decreasing frequencies and during the blocking step one integrates out the high- 
frequency modes above the cut-off fc, keeping the generating functional invariant. We decompose the field variables 
into the high-frequency <j)' x — S\ p \>k^P e%PX ana - low-frequency <Px = Jo<|p|<fc §v e%vx modes. The high-frequency 
modes are integrated out step by step in infinitesimal momentum shells of thickness Afc, 

e -s k -± k {4>] = J x>[<l)']e- Sh ^ +<l> '\ (21) 

with Sk[4>] the blocked action in Euclidean spacetime. The higher- frequency Fourier modes are split further into the 
sum of the saddle point field configuration sp and the remaining field fluctuations: 4>' = cf) sp + ip. In order to evaluate 
the blocked action in Eq. (f2"Tj) , we expand the blocked action in Taylor series around its saddle point field configuration 
and then, for trivial saddle points 4> sp = one arrives at the WH-RG equation [28| 

(2 + kd k )V k = -ilog((l + F fe u )(l + Ff)-(F fe 12 ) 2 ), (22) 

for the LSG model with N = 2 and the second derivatives V k 13 = d^d^.Vk (all dimensionless quantities are denoted 
by a tilde superscript), which is introduced as the sum of the dimensionless mass term and the dimensionless periodic 
potential, 



Vkifa, <h) = M(4>i - fc) 2 + Uh{<k,<h), 



(23) 



7 



which is the generalized Euclidean form of the Hamiltonian in Eq. ([6]) . The ansatz for the periodic part of the potential 
is 

U k =2_j [u«in 2 cos(ni/30i)cos(n 2 /30 2 ) + w„ lM2 sin(rii/30i) sin(n 2 /3(/)2)] • (24) 

ni,H2 

The scale-dependence is entirely encoded in the dimensionless couplings u ni n 2 — u n 1 n 2 /k 2 and v niri2 — v nin2 /k 2 . 
Inserting Eq. (j23|) into the functional evolution equation in Eq. (j22|) and Fourier expanding it we obtain a set of 
coupled differential equations for these couplings. The right hand side of the WH-RG equation in Eq. (|2"2"]l turns out 
to be periodic, while the left hand side contains periodic and non-periodic parts, as well. Separating them, we obtain 
a trivial evolution, 

Mk = Ma , (25) 

hence the dimensionful bosonic mass M k remains constant during the RG procedure. We note that the RG flow 
equation in Eq. ([22]) keeps the periodicity of the periodic piece U k of the blocked potential in both directions of the 
internal space with unaltered length of periods, therefore (3 does not evolve in LPA. 

During the successive integrations of the higher- momentum modes of Eq. (|22|) by the method of steepest descent, 
it is assumed that the saddle points (4>i P , </>2 P ) are zero. It remains trivial till the second functional derivative of the 
blocked action S' k ' (the argument of the logarithm in Eq. (|22[)) is positive definite. When 

S'u' = (1 + V k n ){l + V 22 ) - (VI 12 ) 2 ) = 0, (26) 

the restoring force for the field fluctuations with momenta in the momentum shell k — Ak < \p\ < k vanish, i.e. their 
amplitudes can grow to finite values, and <j>* p becomes non-vanishing, and the spinodal instability appears. Then, the 
tree level blocking relation is needed [1, [U [H, [32| 

Sk-Ak[<f>i,<h]= min (Sk[4>i + 07,02 + <i>2 3 ]) ■ (27) 

<f>l,<f>2 

In case of two layers the exchange of the field variables is a symmetry of the model, therefore the saddle point can be 
considered the same for both fields, namely: 

0f = $ sp = (f> sp = 2p k cos(fcx 1 ) . (28) 

It is assumed that the periodic coordinate dependence appears only in the space direction, which suits well to the 
treatment of the inhomogeneous ground states. Restricting ourselves to a finite interval [-L; L] there is a periodic 
boundary condition for the saddle point, 4> sp {x — L) = <fi sp (x + L). Then the tree level blocking for the boson mass 
M k has the form 

Af fe -Afc(0o - 4>l? = M k {4>l - 2 ) 2 kd k M k = 0. (29) 

so Mfe = M is a constant in the case of tree level blocking as well. Thus, M k ~ fc~ 2 is a relevant coupling for all 
scales. The blocking step 



V k -Ak[0l0° 2 } 



mm 



2 Pk k 2 + 



(u ni n 2 + v nin2 ) cos(/3(ni0? - n 2 0°))( 2 ("i ~ n 2 )Pp k ) 



Kin, - f ni „ 2 ) cos(/3(ni</>? + n 2 02))(2(m + n 2 )(3p k ) 



(30) 



determines the tree level evolution for the couplings. The SG and the MSG models showed that the phases of these 
models can be distinguished by the appearance of the spinodal instability @, [2l[ . Since only the fundamental mode is 
relevant both in the SG and the MSG models, perhaps it is not so surprising that the condition of spinodal instability 
gave a very good approximation for a single coupling. Close to the scale of spinodal instability fcsi the effective 
potential starts to form a parabolic shape [2l|, so the minimum of the blocked potential is situated always at the 
values of the field variables (p\ — 4> 2 = 0. Taking into account the first few couplings (uoij and fin) and inserting 
it into Eq. (|2l)|) one obtains 



(! + 2M k -iioif3 2 -uiif3 2 -iiiP 2 ){l-uoi{3 2 -u 11 (3 2 +i> 11 f3 2 ) = 0. 



(31) 
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The fundamental mode follows the scaling relation 

•fc 2 



"oi = "oi (A) 



A 2 



2M\ 
-2MJ 



/3 2 /16tt 



-2+/3 2 /8tt 



(32) 



The scale fcsi can be situated below or above the mass scale M, so one can distinguish two cases: 

1. kgi < M. The results of Appendix lAl show fin = Un (the massive modes) and Eq. (|3"Tj) reduces to 

(1 + 2M k - uoi/3 2 - 2u 11 p 2 ){\ - u 01 f3 2 ) = 0. 



(33) 



The expression cannot be zero if (3 2 > /3 2 since uoi is irrelevant. If 1 < /3 2 then the expression 1 — iioiP 2 always 
go below zero, because now {t i is relevant and grows up. 

2. ksi > M. Then the massive modes are irrelevant keeping their UV scaling according to Eq. (|32|) but the coupling 
woi scales relevantly when ft 2 < (3 2 , which implies that the expression in the second parenthesis in Eq. (13"T]) can 
be negative. When 1 > /3 2 , then there is no spinodal instability. 

The discussion gives that the critical value of /3 changes as compared to the case of the SG model, now 1 = 167T. 
Since in our model (3 2 = Air the spinodal instability always occurs. It was shown [33| that in this phase the couplings 
u ni n 2 are relevant at least in the UV region. 

We have taken into account 8 couplings when we solved the evolution equations numerically. As in 0, l2lT | we 
obtained that the higher modes do not affect the scaling of the fundamental mode. However the other couplings flow 
by different scaling behavior as obtained by an extended UV RG approach [T3, [33[ . We started the evolution with 
the WH-RG method then at the scale ksi we turned to the tree level blocking equations. In the region (3 2 > (3 2 a 
very interesting scaling law appears, as shown in Appendix [XJ Furthermore, when 1 = Air we run into the region 
of spinodal instability. The numerical results can be seen in Fig.[5J We numerically obtain that the dimensionless 
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FIG. 

M 2 = 



0.01 



5: The scaling of the fundamental mode M01 for two different values bosonic mass, Mi — 0.141 (dashed lines) and 
= 0.04 (solid lines). For Mi the UV (fc -1 ) and the IR (fe~ 3/2 ) scaling laws are drawn as reference. 



fundamental mode goes to a constant as in 2I( independently of the relation between k$i and M , and the IR value 
of fioi(0) is superuniversal, independent of any UV parameter. When ksi < M the flow of the fundamental mode 
qualitatively changes as the scale k goes below M. From Eq. (|32p one obtains Uoi(k — > 00) ~ fc -1 in the UV limit 
but well below the mass scale M when k 2 <C 2M 2 the flow scales according to uoi(k <C M) ~ fc~ 3 / 2 , and after we 
reach the scale ksi and go towards the the IR regime the flow runs into Moi(0). The evolution is similar in the case of 
ksi > M, but since the tree level blocking relation in Eq. (|30|) is independent of M, then below M the scaling cannot 
change. Due to the necessity of the double Fourier expansion we are unable to get as reliable numerical data as in 
[2l| , and we cannot read off the IR form of the effective potential. From the solution of the condition S'£ = 0, [2lT | 
and the qualitative value of Moi(0) ~ 0.14 one expects 

^o = ~(0? + $). (34) 
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The IR physics of the zero density LSG model is free, the dimensionful coupling vanishes as k — ► at 1 = Air. However 
there are two phases with the critical value f3 2 and they can be distinguished by the different scaling behavior of the 
dimensionful couplings, namely when 1 < /3 2 , then the coupling scales as uqi ~ k 2 , while in the case of 1 > [3 2 
the scaling relation is uqx ~ k 13 / 87r . Nevertheless according to the flow of the dimensionless couplings, in the phase 
1 < f3 2 the fundamental mode goes to a constant value, so it is marginal, but in the other phase it vanishes, therefore 
it is irrelevant in the IR limit. 



1. Massless fermions in multi-flavor QED2 



One can easily generalize the previous analysis for multi-flavor systems. We note that the critical value is 1 = 
8nN/(N — 1) '17|, so the spinodal instability always appears at 1 — 4tt. The results of the RG analysis can be easily 
generalized and can show that the dimensionless coupling uj. tends to constant values in the IR limit. Transforming 
the couplings of the scalar model (tik and Mk) into the fermionic ones (gk and nik) one has 

2 1 

M£ = u k = —n 2N m » e«. (35) 

7T 2 

The dimensionful minimal coupling ek does not scale independently of N, therefore the fermionic mass mk should 
follow the irrelevant scaling of Uk 

rn fc ^ ~ fc 2 ^r 0, (36) 

meaning that in the IR limit we always have massless multi- flavor QED2. The pure QED2 with N — 1 is not massless 
since, the tree level blocking gives trivial (which is relevant in d = 2) scaling for the dimensionless coupling Uk so the 
fermionic mass becomes constant, showing that the cases N = 1 and N ^ 1 significantly differs [TH. \u\. 



A. Low density phase 

Now let us turn into the case when b =/= 0. The renormalization is rather involved in finite-density systems due to 
the coordinate dependent ground state field configurations. However in the low density phase average baryon density 
is zero (the low external baryon density is totally screened by the induced density), and the resulting Hamiltonian 
becomes the same as in the zero density case. Therefore we look for the RG evolution of the model eventually among 
constant field configurations just in the case of b = 0, so all the results obtained there are valid in this phase too. 
It implies that the low density LSG model is free, the fundamental mode is marginal, and the dimensionful coupling 
scales as uqi ~ k 2 . The corresponding low density 2-flavor QED2 is massless in the IR limit. 



B. High density phase 



We consider the RG evolution only for the high density limit. Here we always have a non-trivial saddle point cj) sp , 
and the quantum fluctuations takes place around this field configuration. However during the blocking steps <p sp also 
changes, which contributes to the evolution of order O(H ), therefore we concentrate on the evolution only of the 4> sp . 
Using the tree level blocking relation in Eq. (|27|) , one obtains 



Sk-Ak[4>] 



u'k-Ak^ J ^ cos(V2~^) sin(V2^6L) 



2 (di^ sp ) 2 + < cos(%/2^(0 + <j) sp + bx)) 



(37) 



with <j) a constant field configuration. In high density limit the amplitude A\ is small and can be considered pertur- 
batively. Furthermore at high densities one can retain only the fundamental mode. Then the numerically determined 
form of the saddle point field configuration is 

(38) 



cp sp = -A x sin(V27r(&a; + (f)) + Ax sin(V27r0), 
with Ax — u' k /^/2Trb 2 . After identifying the corresponding functionals in Eq. (|38| one obtains 



l k-Ak 



«i 2 



4& 4 



(39) 
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for the infinitesimal blocking step. The blocking relation in Eq. (|39|) clearly shows that dimensionful coupling u' k 
decreases, and tends to zero so as the amplitude A\. Since the ratio Ak/k is kept small in the WH-RG procedure, 
the blocking relation in Eq. (|39p is valid for the dimensionless couplings too, implying that the fundamental mode is 
irrelevant in the high density phase. The quantum fluctuations really wash out the wavy field configurations although 
at that price that in the 1R limit one obtains a free theory in the high density phase too, so the high density 2-flavor 
QED2 is also massless. 



V. CONCLUSIONS 



We determined the phase structure of the bosonized version of the 2-flavor QED2 and showed, that in the tree level 
the induced baryon density totally screens the low external baryon density b. In the case of high densities the screening 
is only partial, furthermore a wavy induced baryon density is obtained. We also performed the renormalization of 
bosonized model in both phases. When b is zero then during the evolution a non-trivial saddle point appears and the 
tree level blocking gives a dimensionful periodic potential which flattens out, implying massless 2-flavor QED 2 in the 
IR limit. It is also the case when b is small, since there is a total screening and the model becomes similar to that of 
zero density one. The evolution of the bosonized model in the high density phase shows that the amplitude A of the 
periodic field configuration decreases as the quantum fluctuations are integrated out step by step so as the fermionic 
mass, and in the IR limit we obtain a massless theory with a trivial ground state. The scaling of the dimensionless 
fundamental mode is marginal in the low density and irrelevant in the high density phases in the bosonized model, 
therefore the phases found at the tree level, survive the RG evolution and exist in the IR limit. 
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APPENDIX A: IRRELEVANT RENORMALIZABLE OPERATOR IN THE LSG MODEL 

The evolution equation will generate the higher order Fourier modes, so one should take them into account according 
to the ansatz in Eq. (|24|) . We introduce a truncation in the number of couplings, namely n\,ni < 2. Then we have 8 
apparently independent Fourier modes. The dimensionless WH-RG equation is 

(2 + kd k )U k = -a log [(l + 2M fe ) + (l + M fc ) (t/ fc n + Uf ) + U^XJf + 2Ul 2 M k - (C/ fe 12 ) 2 ] (Al) 

Now one has to Fourier expand the equation to get the evolution equation for the couplings. The resulting system 
of ordinary differential equation is tackled by a numerical program. As we mentioned the RG-flow always runs into 
the spinodal instability when 1 < /3 2 . In Fig. one can see the typical IR flow of the couplings. One can get the 
following conclusions from the numerical analysis: 

1. We numerically obtained that instead of 8 couplings we have only 5 independent one, Un(k — * 0) = v\\{k — > 0), 
Ux%{k — > 0) = vi2(k — > 0) and U22{k — * 0) = V22(k — > 0), implying that the ansatz of the periodic part of the 
potential in Eq. (|2"5)l can be simply taken as 

_^ 00 
V k = -M (0i - (t> 2 f + cos(tii/30i - naPfa). (A2) 



2 

This result is independent of the value of (3 



n 1 ,n 2 =0 



2. At the scale around the mass of the theory ~ V M k the scaling of the couplings change. There are modes which 
scales trivially below \J M k , namely 

U ni ni ~ k , (A3) 

and they are referred as massive modes. Furthermore the other modes has qualitatively new scaling behavior. 
Numerical results show that in general the Fourier amplitudes scales according to the law 



* 11 in 2 



f,\ni-n 2 \(§^-2)-2S ni , n2 ^ 



11 



>> 

bo 
o 



B 


"be 
o 



-20 



-40 



-60 



-80 



-100 



-120 



-140 











U 1 




- 


■~-^_Uil=V u 












Ul2=Vl2__^^ 











-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 

log(k/A) 

FIG. 6: For /3 2 = 24vr the IR behavior of various Fourier amplitudes is presented for the LSG model. 



with the Kronecker delta 5. 
and El 



, ,,. This result is also valid for arbitrary value of (3, but when < /3 C the Eqs. (|A3[) 
are correct until k > ksj. Below the scale of the spinodal instability the scaling relations of the flow 



changes. 



Starting from different UV initial values of the Fourier modes one obtains that only the fundamental mode £t i 
is sensitive for the initial conditions. The sensitivity matrix has the same structure as in the case of the SG and 
MSG models [|| [2l| where the matrix has only one nonzero column. We numerically obtained that the quantity 



~ni+?i2 

%1 



(A) 



,(*) 



(A5) 



is independent of u„ iri2 (A) and u„ lTl2 (A). Consequently, according to the sensitivity matrix the only relevant 
operator is the fundamental mode uqi which drives the IR behavior of all the other couplings. On the other 
hand according to the scaling law in Eq. (|A4[) fioi ~ k 13 / 87r ~ 2 with positive exponent in the case of (3 > (3 C giving 
an irrelevant scaling. 

By the numerical solution of the RG equation derived for the LSG model we show that for (3 > (3 C one can parameterize 
the RG flow of the couplings by the initial value of the fundamental mode iioi(A), including the IR relevant massive 
modes. The results of the sensitivity matrix shows that the fundamental mode is the only relevant operator of the 
model, nevertheless it goes to zero, so according to common classification of the couplings it is irrelevant. Therefore, 
the scaling of all the higher harmonics and hence the low energy effective theory of the LSG model is driven by an 
irrelevant coupling, showing the essential importance of the investigation of irrelevant operators and the RG technique. 
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